Exact Solutions to Special High Dimensional 0{n) Models, Dimensional Reductions, 
gauge redundancy, and special Frustrated Spin and Orbital models 
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This work addresses models (e.g. potential models of directed orbital systems- the manganates) in 
which an effective reduction dimensionality occurs as a result of a new symmetry which is interme- 
diate between that of global and local gauge symmetry. This path towards dimensional reduction is 
examined in simple 0(n) spin models and lattice gauge theories. A high temperature expansion is 
employed to map special anisotropic high dimensional models into lower dimensional variants. We 
show that it is possible to have an effective reduction in the dimension without the need of compact- 
ifying some dimensions. These models are frustrated and display a symmetry intermediate between 
local and global gauge symmetries. Some solutions are presented. Our dimensional reductions are a 
generlization of the trivial dimensional reduction that occur in pure two dimensional gauge theories. 
It will be further seen that the absence of a "phase interference" effect plays an important role 
in high dimensional problems. By identifying another ( "permutational") symmetry present in the 
large n limit, we will further show how to generally map global high dimensional spin systems onto 
a one dimensional chain and discuss implications. 



I. INTRODUCTION 

In this article, a high temperature expansion is em- 
ployed to examine mappings between problems of various 
dimensionality. By these reductions we will exactly solve 
an anisotropic two dimensional XY and some other mod- 
els. It will be seen that the absence of a "phase interfer- 
ence" effect plays an important role in high dimensional 
problems; this observation will allow us to examine var- 
ious high dimensional models by examining newly con- 
structed one dimensional systems. We will show how an 
effective reduction dimensionality might occur as a result 
of an enhanced symmetry which lies midway between a 
full blown local symmetry and the much more restricted 
global symmetry. Such a symmetry might lead to an 
effective reduction in the dimension without the need of 
actually compactifying any dimensions. By further iden- 
tifying a permutational symmetry present in the large n 
limit, we will further show how to generally map global 
high dimensional spin systems onto a one dimensional 
chain and discuss implications. 

In section ||, We are going to start off by discussing a 
momentum space high temperature expansion. Next, by 
analytically continuing the momenta and making them 
complex (section [II), quaternion (section IV), and of a 
general matrix form we will be able to exactly solve a va- 
riety of two, three, and four dimensional Ising and 0{n) 
models. Many of these models are artificial and corre- 
spond to non hermitian anisotropic long range interac- 
tions. Some of the resulting models will be short ranged 
and of real hermitian form. We will exactly solve some 
special simple two dimensional anisotropic XY models. 
These models may serve as simple caricatures for directed 
orbital models for transition metal oxides (e.g. the man- 



ganates). These models possess a symmetry which is 
intermediate between that of global and local theories. 
Similar models also occur naturally in the low energy 
regime of frustrated pyrochlore and checkerboard anti- 
ferromagnets. 

We will then extend these techniques to the quantum 
case and discuss "2-|-l dimensional Bethe Ansatz" solu- 



tions in section(VI) where we will also solve a special two 
dimensional anisotropic XY model. 



Next, we show in Section(VIl) that the trivial dimen- 
sional reduction of a 1+1 dimensional gauge theory to a 
one dimensional spin chain can be thought of as a trivial 
example of our general mode for dimensional reduct ions. 

After solving these models, in subsection( VIII A ) we 
will apply the high temperature expansion to reduce 
the three dimensional nearest neighbor Ising ferromag- 
net into an illusory formal solvable expression written in 
a determinant notation. In this method, the high tem- 
perature coefficient to each order is arrived at by looking 
a single linear relation originating from a matrix of (su- 
per exponentiall y) incr easing order. 

In subsection( VIII B ) , we will map the nearest neigh- 
bor d dimensional Ising models to a one dimensional 
model whose spectrum is the sum of d mutually incom- 
mensurate tight binding like spectra. This in turn will 
enable us to map the d dimensional Ising model onto a a 
single spin 1/2 quantum particle in 0-1-1 dimensions. 

In section (|x[), we will discuss a "permutational" sym- 
metry present in the 0{n oo) models but violated to 
O (/?'*) in finite n problem. This will allow us to map any 
large n problem onto a one dimensional lattice system. 

One of the main motifs of this article is the absence or 
presence of "momentum interference" effects in high and 
low dimensional problems respectively. Another recur- 
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ring motif is that of a novel symmetry which interpolates 
between global and local symmetry and which effectively 
lowers the dimensionality of the problem. 



and analogously for the 0{n > 1) model 



II. A MOMENTUM SPACE HIGH 
TEMPERATURE EXPANSION 

We will for the most part examine simple classical spin 
models of the type 



(1) 



x,y 



Here, the sites x and y lie on a (generally hypercubic) 
lattice of size N. The spins {S{x)} are normalized and 
have n components. The kernel V{x, y) is translation 
invariant. The n = 1 case simply corresponds to the 
Ising model. 

Later on, when we will discuss quaternion models we 
will replace the scalar product by a slightly more com- 
plicated product. 

The Hamiltonian in the (non-symmetrical) Fourier 
basis is diagonal (/(fc) = Yl,x^i^)^~^^'^'^ ^(^) = 

jfJ2kfi^)^^''^) ^^'^ reads 



(2) 



where v{k) and s{k) are the Fourier transforms of V{x) 
and S{x) respectively. 

For simplicity, we will set the lattice constant to unity- 
i.e. on a hypercubic lattice (of side L) with periodic 
boundary conditions (which we will assume throughout) 
the wave-vector components fc; — ^Ip. where r; is an 
integer (and the real space coordinates xi are integers). 

In the up and coming V = PV. For an invertible 
V{x,y), and an inverse temperature < /? < cxd (via 
a Hubbard-Stratonovich transformation Q , see also Ap- 
pendix A), the partition function: 

Z = Tr exp{-i ^ V{x, y}S{x)S{y)} 



2V, 



= \/det(~^) / l[drj{x)exp[~H^^^'Mx)}] 



(3) 



where 



-El-[(^)"/^-^4/2-i(|ry(f)|)] (5) 

with In/2-i{^) ^ Bessel function. 

The kernel within the dual Hamiltonian H^.'^"-^ is the 
inverse of the kernel V appearing in the original Hamil- 
tonian. It is Coulomb like for the nearest neighbor 
kernel V{x,y) = 

The cosh 77 and Bessel function terms appearing in 
Eqs.(^ and (||) may be viewed as terms of constraint 
securing the normalization of S{x) at every site x. The 
rotational invariance of the scalar product in Eqn.(|l|) is 
manifested by the rotational invariance (|rf| dependence) 
of the Bessel function terms. 

If a magnetic field h{x) were applied, the argument of 
the cosh would be replaced cosh(ry + h) with a similar oc- 
currence (|rf| — > 1 7/4-/1 1) for the Bessel function appearing 
for the 0{n > 1) models. 

The state Potts model can be viewed as a spin model 
in which the q possible polarizations of the spins lie at 
the vertices of a (g — 1) dimensional tetrahedron- in this 
manner the scalar product amongst any two non iden- 
tical spins is cos^^(— l/(g — 1)). Employing this repre- 
sentation, the sum ^-lncosh77(x) may be replaced by 
1^ ^i=i 6xp[S'i • f]\ where the sum is over the q polar- 
izations of the {q — 1) dimensional spin Si. 

For a q state clock model, the Bessel functions are 
similarly replaced by l'^S?=i exp[Si • ff\. For the 4- 
state clock model which we will consider later on this 
is X^ ln[2(cosh77i + cosh 772)]. We will introduce, in 
section (pv|), in the context of the directed electronic 
orbitals appearing in the manganates, "directed clock- 
models" . These models share the same generating func- 
tion J2x^^^i=i ^'^P[^i ' v\ albeit with a polarization (a) 
dependent kernel Vajs {x, y) . 

The partition function of the Ising spins reads 



Z = ^det(-^) j Dije-^" 
det(-^) j X{dr^{x) 

X 

i^77(f)F"i(f, 77)77(7/) 
JJcoshr/(x). 



(6) 



ln[cosh7;(g)] (4) 



For a translationally invariant interaction V{x, y) = 
V{x — y), this may be written in momentum space 
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where v{k) is the Fourier transform 

v{k) = ^ T/(f - y)e*^ = (3v{k). 



(7) 



(8) 



where the sum is over ah lattice sites x. Thus the parti- 
tion function is , triviaUy, 

^ 1 

E (n(^^^"'(-^o)o (9) 

where ( )o denotes an average with respect to the un- 
perturbed Gaussian weight 



exp 



1 
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Y^V-\k)T^{kU-k) 



(10) 



and A/" is a normahzation constant. 

Each contraction {ri{x)rj{y)) (or {ri{k)r]{—k))) leads to 
a factor of V{x, y) = PV{x, y) (or to j3v{k)/N in momen- 
tum space). Thus the resultant series is an expansion in 
the inverse temperature (3. In the up and coming we will 
focus attention on the momentum space formulation of 
this series. All the momentum space algebra presented 
above is only a slight modification to the well known high 
temperature expansions usually generated by the Hub- 
bard Stratonovich transformation directly applied to the 
real space representation of the fields 77(2;) (|J^. It is 
due to the naivete of the author that such a redundant 
momentum space formulation was rederived in the first 
place. However, as we will see, in momentum space some 
properties of the series become much more transparent. 
To make V invertible and the expansion convergent, we 
will shift v{k) by a constant 



v{k) v{k) + A 



(11) 



with A = const > 2d such that v{k) is strictly positive. 
In real space such a constant shift amounts to a trivial 
shift in the on site interaction (or chemical potential) 



V{x, y) V{x, y) + A.5s,y 
For asymptotically large |?7| the nontrivial 
^duai = _^in|coshr/(f)| 

X 

is linear in |7;|. The Gaussian generating 



(12) 



(13) 



fjdual * 



2 E'^(^)^ ^(^^y)'niy) 



(14) 



is a positive quadratic form and dominates at large rj. 
Note that the convergence of this series in Eqn.(^ is 
better than that appearing in the canonical field theories 
where the perturbing Hi dominates over Hq for large 
fields and, as a consequence, all canonical expansions are 
blatantly asymptotic [||. 

So far all of this has been very general. Now let us con- 
sider the d dimensional nearest neighbor Ising model. To 
avoid carrying minus signs around, we will consider an 
Ising antiferromagnet whose partition function is, clas- 
sically, identically equal to that of a ferromagnet. The 
exchange constant J will be set to unity. The correspond- 
ing momentum space kernel 



v{k) ^ A + 
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1=1 



,iki 



-ik. 



(15) 



The resulting high temperature series (In Z = 
~NJ2p fpP^) must be identical to the one derived by 
conventional real space methods which is convergent for 
all (3 < Pc (the inverse critical temperature): we retiter- 
ate that this convergence proves that the series we obtain 
is not asymptotic and free from many of the pathologies 
of common diagrammatic expansions 

As the reader might guess, the uniform shift A only 
leads to a trivial change in Z (i.e. to a multiplication 
by the constant exp[— /^AA^]). In the forthcoming we will 
set A = Q. 

Owing to the relation 



dk e 



IkN' 



2ttS 



N'.O 



(16) 



all loop integrals are trivial. 

All stated here and below also holds exactly for finite 
size (of sides L > 1) systems with periodic boundary 
conditions. Here 



^exp[ifc„7V'] = LSn',0 



(17) 



replaces Eqn.(p^). In Eqn.(p7|) the sum is over all kn = 
with —7T<kn < TT. Notc, howcvcr, that Eqn.([l7|) is 
far more general and holds for a large class of finite size 
systems. [0. 

In the thermodynamic limit, a simple integral is of the 
type 



M 



E 



(18) 
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(when n = 1 this integral is 2d etc). A related integral 
reads 



(27r)" 



E 



M 



(19) 



spatial components I satisfying the delta function con- 
straints. In a given diagram not all propagators need to 
correspond to the same spatial component fc' . The com- 
binatorics of "in how many ways loops may be chosen to 
satisfy the latter delta function constraints?" boils down 
to the standard counting of closed loops in real space. 
Explicitly, in general d, 



'j<n i- "H<M — —ITT-i^M 6tC. 



where in the last sum m'^^ ^ «_j„ 
Note that the high temperature expansion becomes sim 
pier if instead of symmetrizing the interaction (each bond 
being counted twice- once by each of the two interacting 
spins) one considers 



v{k) — 2^^exp[iA:i]. 



(20) 



1=1 



Here each spin interacts with only d (and not 2d) of 
its neighbors separated from it by one positive distance 
along the d axes ef^j^. The factor of 2 originates as each 
bond is now counted only once and therefore the corre- 
sponding bond strength is doubled. The partition func- 
tion Z is trivially unchanged. 

In the most general diagram the propagator momenta 
Wa} (1 < a < number of propagators) are linear combi- 



loops 



nations of the independent loop momenta {fcbjj^]^ 
qa = Mabh 



(21) 



with integer coefficients Mat- 

Symmetry factors aside, the value of a given diagram 
reads 

(2/3) 



1 = 1 b 

propagators 



loops 

n 4^E-p[^E^^-fc^]) (22) 

b=l ^ ' 1=1 a.b 



- i.e. upon expansion of the outer sum, just a product 
of individual loop integrals of the type encountered in 
Eqn. (|l^) for each k^. 

The latter integral in d = 1 reads 



loops propagators 

/= n '^( E ^'^-",0) 

b=l a=l 



(23) 



with a Kronecker delta. The situation in d > 1 becomes 
far richer. One may partition the integral into all pos- 
sible subproducts of momenta corresponding to different 



all partitions of {qa} into d sets Ci 



d 

n 

1=1 



loops 

n S{J2Mab,0). (24) 

6=1 aeCi 

Note that entirely identical expressions would be 
reached for 0{n > 1) models- the diagrams to be drawn 
and their integral expressions s are exactly the same as 
for the Ising system. The sole difference is that in the 
0{n > 1) case, the vertices of various orders (given by 
the Taylor expansion of the Bessel functions) will have 
different weights. 

We may imagine coloring, for each individual term in 
Eqn. (p^) , the different sets Ci of propagator lines {g^} 
corresponding to different dimensions by different colors 
to produce a graph colored with d colors. It is amusing 
to note that, by the four-color theorem, the most general 
topology of planar diagrams is present at d = 4 dimen- 
sions. 

Some simple technical details of lattice perturbation 
theory vis a vis the standard continuum theories are pre- 
sented in p5[ |. 



III. AN EXACT SOLUTION TO A SPECIAL 
THREE DIMENSIONAL ISING MODEL 

If the spin spin interaction on a cubic lattice has a non 
hermitian kernel of the form 

V{xi,X2,X3) = -^x2,O^X3,o(^a;i4 + Sx^-l) 

1 sinhTi-A r 

(-1) X2 , 2 ^^1:0 A((5^2,l -f (5:^2,-1) 



+ 1X3(6x2,1 - Sx2-l) 



(25) 



then, for all real A, the Hclmholtz free energy per spin 
/(/i = 0,/?>ln(l + \/2)/2) is 



where 



Pf = -ln(2 cosh 2/3) - 
^ d01ni(l + 0~~^2~2^ 

_ 2sinh2/3 
'^^ cosh^ 2/3 



(26) 



(27) 
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and the internal energy per spin 
r 2 

u = -coth2/3 1 + -k'Ki{k) 



(28) 



to write the free energy as the free energy of an 
anisotropic two dimensional nearest neighbor ferromag- 
net with the kernel 



where Ki{k) is a complete elliptic integral of the first 
kind 



V{xi,X2) = -~Sx2,ni^xi,i + Sx^-i 



(36) 



tt/2 



and 



2tanh2(2/3) - 1. 



(29) 



(30) 



i.e. a nearest neighbor Ising ferromagnet having the ratio 
of the exchange constants along the xi and X2 axes equal 
to r — J2/ Ji (throughout we set Ji to imity). Introduc- 
ing 



1 



The proof is quite simple. The expressions presented 
are the corresponding intensive quantities for the two di- 
mensional nearest neighbor Ising model 0] . 

Let us start by writing down the diagrammatic ex- 
pansion (absolutely convergent for T > Tc) for the two 
dimensional Ising model (where Tc = 2/ ln(l + \/2) and 
set 



v{qi,q2) v{qi,q2 + i\qz) 
in all integrals. This effects 



(31) 



sinh 2/3 sinh 2(ir 

and defining 

F{e) = ln{2 [cosh 2/3 cosh 2/3r 
+u-^{l - 2ucos6' + 1*2)1/2]}^ 

we may write the free energy density as 

f3j = -^ r Fie) de. 

27r In 



(37) 



(38) 



(39) 



dfcae-^'^^^"'"' (32) 



for all individual loop integrals. The integrals are nonva- 
nishing only if {Nb = 0} for all loops b. When N'-°°p = 
for a given loop then the second (fca) integration leads to 
a trivial multiplicative constant (by one) . In the original 
summation 



Note that the same result applies for the "ferromag- 
netic" 

V{xi,X2, X3) = -Sx2,oSx3,o{<>xi,l + Sxi-l) 

1 



'JSxi^Q 



E 



(27r)° 



(33) 



X'^ + xj 
+i{6x 



{{3x2,1 + Sx2-l) 

.,-1 - 0x2,1)^ s- 



(40) 



and thus an additional factor of L is introduced for each 
independent connected diagram. Thus the free energy 
per unit volume for a system with the momentum space 
kernel v{ki,k2 + iXk^) is the same as for the two di- 
mensional nearest neighbor ferromagnet with the kernel 
v{ki, k2). Fourier transforming v{ki, k2 -f iA/ca) and sym- 
metrizing [V{x) + V{—x)]/2 V{x) one finds the com- 
plex real space kernel V{xi,X2,X3) presented above |^ 

For the more general 

V{xi,X2,X3) = ~Sx2,oSx3,o{Sxi,l + Sxi-l) 

(_1).3 



The proof is simple- the partition function is identically 
the same if evaluated with the kernel in Eqn. (^) instead 
of that in Eqn.(M) if the spins are flipped 



Six) (-l)^-'S'(f). 



(41) 



What about correlation functions? The equivalence of 
the partition functions Z[P, h{xi, X2)] implies that in the 
planar directions the correlations are the same as for two 
dimensional ferromagnet. 

Changing the single momentum coordinate 



-J^Sxi.O 



A2-|-a;2 



{{Sx2,l + 6x2-1) 



effects 



we may define 



-i{6x2~i - ^2:2,1)^} 



ttJ 



2A sinhvrA 



(34) 



(35) 



v{ki,k2) v{ki,k2 + iA/ca) 



G{ki,k2) G{ki,k2+iXk3) 



(42) 



in the momentum space correlation function. The alge- 
bra is identically the same. 
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The connected spatial correlations (for T ^ Tc) are 
exponentially damped along the X2 axis as they are in 
the usual two dimensional nearest neighbor ferromagnet 
(see for the two dimensional correlation function). 

Of course, all of this can also be extended to one- 
dimensional like spin systems. 

The two dimensional spin-spin kernel 



\ 1/ sinhTrA 



ix2{5x2A - Sx2~i)] (43) 



leads to one dimensional behavior with oscillations along 
the X2 axis. More precisely, the momentum space corre- 
lator for any 0{n) system in one dimension reads 



G(fc) 



+ 



1 



,ik+r _ Y ' 



(44) 



For the Ising {n — 1) system r = — [lntanh/3]. One 
may set k = ki + iXk2 and compute the inverse Fourier 
transform. The correlations along X2 are oscillatory with 
temperature dependent wavevectors. In other words, the 
effective correlation length along X2 is infinite. 

If the two dimensional spin-spin interactions in 
Eqn.(|4^) are augmented by an additional on-site mag- 
netic field h then the two dimensional partition function 
reads 



Z = Ai' + A 



where, as usual. 



A± = cosh Ph ± \Je^P sinh^ ^/i - 



-2/3 



(45) 



(46) 



That "complexifying" the coordinates should not 
change the physics is intuitively obvious: if one shifts 
fci — > fci + Afc2 = k'l with real A in the one dimensional 
v{ki) then the resulting kernel v{k'i) plainly describes a 
stack on chains parallel to (1, A); the spins interact along 
the chain direction yet the chains do not interact amongst 
themselves. The free energy density should be identically 
equal to that of a one dimensional system- no possible de- 
pendence on A can occur. All that we have done in the 
above is allow A to become complex. 

In 



v{ki,k2) — > v{ki,k2 + Afcs) 



(47) 



the correlation functions along a direction orthogonal to 
the direction (1,A) in the (x2,a;3) plane vanish. If A is 
complex then there are no correlations along an orthog- 
onal direction in the "complex" space. 

Formally, all this stems from the trivial fact that the 
measure dkidk2 = const(dkdk*), whereas v{k) depends 
only on the single complex coordinate k. 



Other trivial generalizations are possible. For example, 
we may make both fci and ^2 complex in v{ki,k2) to 
generate a four dimensional spin kernel 



t(<- 



^^^^ sinhTrAi ^^ 
Ai + xl 



^^ sinh^Aa 



F(a;i, 2:2.2:3, a;4) 

Al (^0:2,1 + Sx2,-l) 
iX3iSx2A - '>X2-l 

A2(^a;4,l + Sxi-l) 



(48) 



which leads to the canonical (nearest neighbor) two 
dimensional behavior for arbitrary Ai and A2. 

All of this need not be restricted to cubic lattice mod- 
els. We may also analytically continue v{k) of the trian- 
gular lattice etc. 

For the triangular antiferromagnet/ferromagnet 



k \/3 

v{k) — ±2[cos ki + cos(— H — ^^'^) 

, ki \/3 , 
+ cos{-— + —k2)]- 



(49) 



Making ^2 complex leads to interactions on a layered 
triangular lattice etc. 



IV. HIGHER 0{N > 1) MODELS AND 
SYMMETRIES INTERMEDIATE BETWEEN 
GLOBAL AND LOCAL 

The multi-component 0{n > 1) spin models display 
a far richer variety of possible higher dimensional exten- 
sions. 

The constant A need not be only a complex number, 
it may also be quaternion or correspond to more general 
set of matrices. 

Let us examine the extension to "quaternion" k (or 
with the matrices ifii, i(T2,and ia^ (with {ffa} the Pauli 
matrices) taking on the role of i,j and k). If we were 
dealing with an 0(2) (or XY) model, the same high tem- 
perature expansion could be reproduced: 

We may envisage a trivial extension to the usual scalar 
product Hamiltonian 



H 



x-y)Six) ■ S{y) 



(50) 



X:y 



to one in which the kernel is no longer diagonal in the 
internal spin indices a, P = 1,2. 



H = ^^Va,i3{x - y)Sa{x)S,3{y)- 



(51) 



6 



The Hubbard Stratonovich transformation for Z = 
Tr{exp[-/3i?]} with the bihnear H = lY.^pSJ^KpSp 
proceeds as before. The super-index n = {x, a) now la- 
bels both the physical coordinate x and the polarization 
a. And, just as before, the trace of exp[iSKTiK.] over all 
{Sk} is rix ^n/2-i(l^(2?|)- The high temperature expan- 
sion of section may be reproduced word for word with 
the kernel Va^f3{x — y) now generated by each contraction 
of 770,(2;) and rii3{y), or in momentum space Va,f3{k)Sj:^j:, q 

is generated by the contraction of 770 (fc) and rjfj{k'). In 
the evaluation of the free energy or the partition function, 
a given bubble diagram reads 



loops 

n 

6=1 



^rr{exp[z^Mafcfc[]} 



(52) 



1=1 



a,b 



where, as before, a spans the number of propagators, 
b spans the number of loops, and fc^ with 1 < Z < d 
denotes the l-th component of the d— dimensional wave- 
vector within the b-th loop kf,. 

Let us start off with a nearest neighbor one dimensional 
XY chain and consider the transformation 

exp[ifci] exp[i(fci -I- iAifc2(Ti + i\2kzCT2 + i^^k^a^)] 

= exp[7fci — w ■ a]. (53) 

In the argument of the exponential the identity matrix 1 
commutes with w ■ a and the exponential may be trivially 
expanded as 

exp[7fci — w ■ (j\ ~ cxp[jfci] exp[— u; • a]. (54) 

Once again in higher dimensions (this time d = A) the ki 
integration will reproduce the familiar 



dki exp[iiV'fci] = 2t:5n',o 



(55) 



for each independent loop momentum k and the integra- 
tion over the remaining fc2,^3, and k^ components will 
simply yield multiplicative constants. 
Thus the kernel 



v{k) — exp[ifci — w ■ a] + h.c. 



(56) 



will generate a four dimensional XY model which has the 
partition function of a nearest neighbor one dimensional 
XY chain. The corresponding real space kernel 



Val3ixi,X2,X3,X4:) 



" d'^k ^ / 

——J exp[ifc • x] ( exp[7(fci 

+i\ik2CFi + iX2kz(T2 + i\^kAa3)]ai3 + h.c^ . 

In general, this kernel is no longer diagonal in the internal 
spin coordinates (unless w happens to be oriented along 
0-3) ■ 



.Xik2 A2fc3 A3A:4 

w = (-p^, -p^, -p^ 

\w\ \w\ \w\ 



and 



exp[— (T • w] = cosh(|?«|) — sinh(|it;|)((T • w) 



(57) 



(58) 



For general {Xi}^^i the inverse Fourier transform is non- 
trivial. The scalar product is replaced by a more compli- 
cated product amongst the components. The real space 
interaction kernel is once again algebraically long ranged 
along the a;2,a:2 and X4 axes. 

Just as quaternions (or Pauli matrices {ffi}) may be 
employed we may also consider the transformation 



fci fci + k^+2l^ 



(59) 



with the Dirac matrices {7''}, with its associated propa- 
gator exp[j(fci + fc^+27'')] to produce solvable models in 
2 < d < b dimensions (depending on how many momen- 
tum components are contracted with the gamma matri- 
ces) . The partition function of these models is identically 
that of a nearest neighbor one dimensional system. An 
infinite number of variations along this theme are possi- 
ble. 



In this manner we may generate a multitude of solv- 
able, slightly more realistic, models with a real hermi- 
tian kernel. 

The hermitian kernel 



v{k) = —2 cos(fci + /e2<T3) 



(60) 



leads to the partition function of a nearest neighbor one 
dimensional 0{2) chain. 

The proof amounts to a reproduction of the calcula- 
tions given above for the four dimensional 0(2) system. 

In real space the latter kernel (Eqn.(^)) leads to the 
0(2) Hamiltonian 



H 



(ij) along (1,1) 
(ij) along (1,-1) 



(61) 



where the n = 1 (or x component) of the spins interact in 
the first term and only the y-components of the XY spins 
appear in the second term; the indices i and j are the two 
dimensional square lattice coordinates. The first term 
corresponds to interactions along the (1, 1) direction in 
the plane and the second term corresponds to interactions 
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along the (1,-1) diagonal. The two spin components sat- 
isfy + [S'f^'']^ = 1- i.e. are normalized at every site 
i. By our mapping, this model trivially has the partition 
function of a one dimensional nearest neighbor 0(2) spin 
chain. Thus its free energy per site / = — ln[/o(/3)] 
where Iq denotes the associated Bessel function. 

It amusing to trivially note by inspection that this ex- 
plicitly works for works the periodic 2x2 system on the 
torus where the two diagonals decouple and along each 
diagonal one has the interactions of the standard 0{2) 
symmetric one dimensional XY model. 

This system is frustrated- glancing at the Hamiltonian 
of Eqn.(|6l|) we note, that it is impossible to saturate the 
bonds associated with both the x and y polarizations 
of the spins simultaneously. The state in which the x 
component interactions are minimized is geometrically 
orthogonal to one in which the y— component interactions 
are minimized. In fact, we just saw that the frustration is 
so large that it thwarts the generic low temperature two 
dimensional quasi (algebraic) long range order and leads 
to a canonical one dimensional like behavior instead. 

We may also generate interactions along arbitrary 
tilted rays by considering 



v{k) ~ — 2cos(/ci + Afc20'3) 
to show that the Hamiltonian 

H = - 

(ij) along (1,A) 
(ij) along (1,-A) 



(62) 



is equivalent to 



H - -^Si ■ Sj. 



(63) 



(64) 



By tuning A ^ 0, the interactions along the two "clap- 
ping" directions (1,±A) fold back and degenerate into 
the original one dimensional, 0(2) symmetric, chain. 
Similarly, by choosing the diagonal 



Here as all v matrices are diagonal in the internal 
spin indices, the traces of the two dimensional system 
(Eqn.(^)) are identically the same as those of the one 
dimensional system of Eqn. (|67|) . Schematically, in the 
two dimensional system, for each given loop momenta 
having the two spatial components (fcj, k^^), the integrals 
are exactly the same as they are for the one dimensional 
case: 



Tr 



dk\ 



dk\ 
"2^ 



„ik^N' 





AklN' 



Tr 



dk\ 
"27 








After integrations, each one of the identical diagonal 
entries will be given by Eqs.(p4|). These are the exactly 
the same integrals that we have evaluated for the Ising 
case. (As told, it is merely the weights of the vertices 
which change in the general 0{n) model.) The dimen- 
sional reduction of Eqn.(|6^) is ambivalent to the special 
nearest neighbor short range (or even lattice character) 
of the problem. If the x-components of the spins interact 
with a kernel which is xi direction dependent and the 
y-components would interact with the same kernel along 
the X2 direction then a dimensional reduction would be 
possible. In general, and in the continuum in particular, 
one may, of course, perform rotations within the exter- 
nal momentum space and internal spin space in order to 
examine dimensional reductions for all interaction direc- 
tions. 

As told, on the lattice, all of this need not apply to 
special short range interactions. More generally, we may 
replace ki ki + fc20'3 within an arbitrary kernel 



R 



v{ki) — —2 2, Vr cos(rfci) 



(68) 



of spatial range R> 1. Employing Eqn.([55|) once again 
we note that any XY spin system with a Hamiltonian 



v{k) = -2cos(fciP+ + k2P-), 



(65) 



with the projection operators P± — i(l ± cts), we may 
prove that 



H=- E sr^s^' 

(ij) along (1,0) 
(ij) along (0,1) 



is equivalent to the one dimensional 
H = - ^Si ■ Sj. 



(66) 



(67) 



(ij) along (1,1) 



E ^^s^'s^' 



(69) 



(ij) along (1,-1) 



with a general translationally invariant Vij = f{\i — j\) of 
an arbitrary range R (i.e. one which does not necessarily 
link only nearest neighbor diagonal sites i and j) has 
exactly the same partition function of a one dimensional 
system whose Hamiltonian 



IJ Si • Sj 



(70) 
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with the arbitrary range interaction kernel Vij . 

Similarly, the partition function of an XY system hav- 
ing the three dimensional kernel v{k) = — 2cosfci — 
2cos(A;2 + fcatTa) is equal to that of the two dimen- 
sional nearest neighbor XY ferromagnet which exhibits 
Kosterlitz-Thouless like behavior at sufficiently large in- 
verse temperature /3 (or temperatures T < Tkt)- Here 
each of two spin components interacts within a differ- 
ent subplane. The coupling between the two spin com- 
ponents due to the normalization constraint apparently 
plays no role in leaving the system two dimensional. 

In a similar manner, any d! dimensional spin system 
(at temperatures T > Tc) can be made have an effective 
dimensionality 1 < d < d' without the physical need of 
actual geometrical compactification. 

The models above display special symmetries, interme- 
diate between those of the standard two spin 0{n) models 
which display a global group symmetry and lattice gauge 
theories which a local group symmetry at every lattice 
site. Note that the transformation 



qixi,X2) qixi,X2) 

^a=l 'lxi'->a=l 

qixi,X2) q{xi,X2) 

^a=2 ~^ '1x2^ a=2 



(71) 



with the (2L) arbitrary, "gauge", degrees of freedom 
rjx^ — ±1 leaves Eqn. (|66|) invariant. Here {xi,X2) de- 
note the two spatial components of the spin coordinate 
X. Note that by this symmetry, any particular state (in- 
cluding the ground state) is, at least, 2^^ degenerate- 
yet another manifestation of "frustration" . The effective 
gauge redundancy of Eqn.(^) is midway between a full 
blown gauge theory (where redundant gauge degrees of 
freedom live on every site) and a globally symmetric the- 
ory where only one global symmetry exists. The number 
of gauge degrees of freedom in Eqn.(|7l|) is the same as 
that in a one dimensional Z2 gauge system. In general, 
the manifold of the gauge degrees of freedom can assume 
any effective dimensionality in between the two standard 
extremes- d (gauge theories) and zero (globally symmet- 
ric spin models). The 0(4) spin model 

H=- J2 (^f^^f^+^f^f) 
(y) along (1,1) 

- E (S^'s^^+Sl'^sf) (72) 
(y) along (1,-1) 

displays (2L) independent 50(2) (or L''(l)) gauge degrees 
of freedom ry residing on each individual diagonal. Ex- 
plicitly 

{q{xi,X2) . ■q(xi,X2) /q{xi,X2) , ■qixi,X2)-s 

qixi,X2) ■ qixi,X2) f q{x 1 ,X 2) ■ q{x i ,X 2) ^ y ) 



with arbitrary phases along the 2L diagonals di^2, 
is a symmetry. By our mapping, this two dimensional 
system with ~ L redundant 0(2) gauge degrees of free- 
dom is equivalent to an isotropic, nearest neighbor, 0(4) 
spin chain having only a single zero dimensional global 
symmetry. The variations are endless. 

Although we have focused attention on equivalence to 
simple models with rotationally invariant scalar products 
(e .g. Eqn.(f70|)), much of what we said is also generally 
valid for non- isotropic actions [ pTf . 

Finally, in order to set the stage for the discussion of 
a cartoon of directed electronic orbitals (e.g. p or d or- 
bitals appearing in some systems (e.g. the manganates), 
we will now introduce "directed clock models". These 
models are just the standard TV— state clock models with 
the additional twist that an interaction ( "hopping" ) can 
occur only along an axis parallel to the spin direction (a 
cartoon for the direction along which the electronic or- 
bital is extended). The high temperature diagrammatic 
expansion for these models may proceed as before, with 
the provision that the generating function (which merely 
sets the weights of all vertices) is no longer a rotation- 
ally invariant Bessel function but rather X^i^i 6xp[Si • ff\. 
For the directed 4-state clock model, the four possible 
spin orientations are (±1,0) and (0,±1) and the Hamil- 
tonian is none other than that of Eqn. ( |66| ) . Our mapping 
will show that this is equivalent to the standard (undi- 
rected) one dimensional nearest neighbor 4-state clock 
model whose partition function is trivially twice of that 
the nearest neighbor one dimensional Ising model. The 
generating function is 

(JJ[2(cosh77i -I- coshr/2)]) 

X 

(n[2cosh(,7+)])(n[2cosh(r;_)]), (74) 

X X 

with ri± = [rji ± r]2)/2. The partition function is a direct 
product of the partition function terms for two indepen- 
dent Ising systems. This has a much clearer intuitive 
meaning (see Appendix B)- the direct product space of 
two Ising spins sitting at every site gives rise to two di- 
mensional spins whose components are (±1,±1). Upon 
rotation and normalization, this becomes the four state 
clock model just discussed. 

In general, it is easy to see that the partition function 
of the (2^*) state directed clock model (in which each spin 
has d components) embedded in an external space of d 
dimensions is equivalent to that of the one dimensional 
Ising model (more precisely, d replicas of Ising systems) . 

We may similarly address ||l^ the "squared directed 
4-state clock model" 



H 



(l)^(l)]2 



(ij) along (1,0) 
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(ij) along (0,1) 

to show that its partition function is equivalent to that 
of a one-dimensional model with 

H = -Y1 ^f^f+i = (= -^)' (76) 

which is trivial. 

The reader is referred to Appendix B for a perspec- 
tive on the conventional, non-directed, and directed clock 
models in n > 2 internal spin dimensions. 

V. POTENTIAL PHYSICAL REALIZATIONS 

Some frustrated systems might have realizations close 
to the special symmetry points at which we proved the 
occurence of exact dimensional reductions. 

A. The Manganates 

Before giving a superficial flavor of one potential sys- 
tem where models and symmetries such as those intro- 
duced in the previous section might be of real practical 
relevance, let us quickly outline the logic. As seen in the 
central figure of Fig.(|l|), orbital wave functions can be 
cigar like and highly direction. In such a case, overlap 
integrals will be much larger for hopping or interactions 
along one axis than along another. Hopping and/or spin- 
orbit couplings might depend on the state of the orbital 
direction ( "polarization" ) . 

Polarization dependent interactions akin to those in 
Eqn.(|6ll) may be of relevance in discussing directed or- 
bital problems in two dimensions (e.g. those associated 
with the electronic orbitals in transition metal oxides of 
perovskite structure) [jl^ . The non-isotropic electronic 
orbitals lead to overlap terms which are strongly direction 
dependent. In these materials (e.g. La2_2:Sr2,Mn03), the 
surrounding oxygen atoms introduce crystal field split- 
ting of the d orbitals of the caged transition metal (e.g. 
Mn). Wave functions pointing towards the oxygen ions 
may attain higher energy than those that point in be- 
tween. The two "bad" higher energy d orbitals are 
jx^ — y^) and |3z^ — r^) (customarily denoted as Cg or- 
bitals) are depicted in Fig.(|l]). The ground state is deter- 
mined by Hund's rule. A Mn'^+ ion, for example, has a d^ 
configuration and consequently, one electron populates 
on the states shown in Fig.(|l]). These states dominate 
the orbital physics. 

The directed four or six state clock models may be re- 
garded as oversimplified cartoons for the orbital lobes in 
systems such as these. The overlap terms allow hopping 
only in the direction lying along the axis of symmetry of 
the orbital. They are a very crude model for 




FIG. 1. Left and middle: the conventional real space func- 
tions (f)a = — y^, — Sz^ — for the cubic E dou- 
blet. Right: the complex combinations (f>a^i(l>b are octupolar 
states. 

H = Y.^^Ma'^,<y + h.c.) (77) 

where the hopping amplitudes {iy } depend on internal 
indices indicating the relative orientation between the 
orbitals at sites i and j. The hopping amplitudes be- 
tween two neighboring Mn atoms depend on the overlap 
of the d-orbitals with the oxygen p-orbitals which lie in 
between. 

Another origin of such a behavior may result from spin- 
orbit couplings. 

B. Frustrations and Transitions within the 
checkerboard and Pyrochlore lattices 

As demonstrated by E. Berg et al. the S = 1/2 
Heisenberg antiferromagnet on the frustrated checker- 
board and pyrochlore lattices, can be related to models 
of the form similar to those that we have looked at till 
now. 

For instance ||l^, in the low energy regime, the Heisen- 
berg model on the checkerboard lattice shown in Fig.(||), 
is related to the Ising like model 

H = -JY,{S, ■ e,,){S, ■ %) -hY^St. (78) 

Here the directors may point along directions at an- 
gles 7r/3(— 7r/3) away from the x axis for horizontal (ver- 
itcal) bonds respectively. The spins correspond to ver- 
tical or horizontal dimers Had the directors 
been along orthogonal axis, we would just obtain a triv- 
ial counterpart of the Hamiltonian of Eqn.(|6l|). At the 
orthogonal point (if the directors were orthogonal), the 
system is exactly be one dimensional. Sans the applied 
mangetic field, the free energy density at this point is 
—(3~^ In cosh/? times a combinatorical factor equal to the 
logartihm of the number of ways in which the plane may 
be tiled with horizontal and vertical lines such that each 
point belongs to exactly one line. 

A deviation of {e^ } from orthogonality, leads to the re- 
moval of the degneracy associated with the symmetries of 
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FIG. 2. The checkerboard lattice. The two singlet ground 
states of the crossed plaquettes may be represented by an 
Ising pseudospin lying along the horizontal axis. The effec- 
tive interaction amogst these pseudo-spins is encapsulated by 
Eqn.©. 

Eqn. (|7l|) and its likes, and the system is indeed expected 
to display canonical two dimensional behavior and have 
a transition at a finite temperature. 

A similar model may be written for the pyrochlore lat- 
tice @. 

Thus, some geomertically frusrated models are close 
to yet removed from the special symmetry points where 
dimensional reduction would occur. 



VI. SPECIAL TWO DIMENSIONAL BETHE 
ANSATZ SOLUTIONS 

We may consider quantum (spin, fcrmionic, bosonic) 
extensions in which the unperturbed action Sq links two 
nearest neighbor sites. By the Trotter formula, the par- 
tition function may be written as a lattice sum over a 
(d + I) dimensional system of thickness f3 (the inverse 
temperature) along the imaginary time axis (the (d+l)- 
th direction). The (d+l)-th momentum component kd+i 
is Matsubara frequency lo whose allowed values {a;„} are 
dictated by the periodic/antiperiodic boundary condi- 
tions along the imaginary time axis that must be im- 
posed for a bosonic/fcrmionic system respectively. The 
Green's function conjugate to Sq in the undualized (non 
Hubbard-Stratonovich transformed) system reads 



Go — 



d-l-l 

1=1 



(79) 



If the magnitude of the squared mass |tooP = 1^1 = 
const > 2{d + I) then a geometric series expansion 
in may be performed to generate a sum of har- 

monics. Within each diagrammatic term, we may take 
ki ki + i\k2 or other similar extensions to create high 
dimensional problems that have the same partition func- 



tion as that of a solvable canonical one dimensional sys- 
tem. 

For a relativistic spin-1/2 system, we may employ the 
naive fermionic propagator 



Go = [/3(7m' 



(80) 



and once again expand the resulting diagram in powers 
of m^^ to obtain a sum of harmonics. This holds for all 
canonical lattice theories in which the kernel are periodic 
functions of the momenta and thus, when well behaved, 
may be expressed as sums of harmonics. 

Eqn.(|55|) forces all diagrams with substitutions such as 
fci ^ fci + i\k2 to attain the same value as they would 
attain in the original system. 

If an exact Bethe Ansatz solution is known for a 
fermionic one dimensional system with nearest neighbor 
hopping then in some cases it may be extended to two 
dimensions where the hopping matrix element K{xi,X2) 
for separation x = {xi, X2) attains exactly the same form 
as the two-dimensional kernel V{xi^X2) just given pre- 
viously. Loosely speaking, if in a fictitious electronic 
system the hopping matrix element would be the non 
hcrmitian 



K{xi,X2) 



1/ sinhTrAi 



Ai((5x2 



i) 



-1-1X3 ((5: 



2:2 ,1 



)]) 



(81) 



then the system would essentially be a one dimensional 
Luttinger liquid along xi . If the kernel would instead de- 
cay algebraically along the imaginary time axis the prob- 
lem would become that of dissipating system. 

A hermitian hopping matrix element could also do the 
trick if we were to consider polarization dependent hop- 
ping (or "interactions") analogous to those in Eqn.(|6l|) 
and slightly more complicated variants. As noted ear- 
lier, such models and the more general equivalence of the 
two dimensional system of Eqn. (|69| ) with the one dimen- 
sional one of Eqn. (^o|) may be of relevance in discussing 
directed orbital problems in two dimensions in systems 
like the manganates |l^. 



VII. DIIVIENSIONAL REDUCTION IN LATTICE 
GAUGE THEORIES 

Similarly, the standard lattice gauge theories may be 
mapped to higher dimensional lattice gauge theories. Let 
us further show also the opposite- the dimensional reduc- 
tion that trivially occurs in going from a two dimensional 
gauge theory to a one dimensional spin chain is one spe- 
cial case in our general way of constructing models of 
different dimensions that are equivalent. 
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A pure lattice gauge theories (for any group) may be 
mapped onto a Coulomb gas model with nontrivial onsite 
interactions. The advantage is that in the Coulomb gas 
all interactions are between pairs and our machinery for 
dimensional reduction can be brought to bear. To avoid 
carrying too many indices around we will now consider 
the Abelian Zi (Ising) group. 

Let us consider the even sublattice sites 



Tr 



O(mod 2) 



(82) 



in a hypercubic lattice, where Xa labels the a-th compo- 
nent of the position vector of a site on the lattice. On 
this sublattice dual fields are introduced. 

Each bond variable C/y is connected to one and exactly 
one lattice site (i or j) which belongs to the even sub- 
lattice (where the fields C exist). The partition function 
may then be schematically written as 



Z = I DUexp[Kj2uUUU] 
J DUY[ expi^Vo^p.^sU^UpU^Us] 



N DU i?Cexp[-lC/^,Ca/3C75 - KcpUo^Up] (83) 



where the kernel V„ 



^'Q/3;75 = K^\ap--iSls/2 where the Greek 
indices denote the bonds and a(3 denotes the site common 
to bonds a and (3. The distance between two sites was 
marked by |a/3 — ^5\. At each sublattice site i there 
are 2d{d — 1) fields {Ca/s}- As each link Uij has only a 
single end point i which lies in the even sublattice, only a 
dependence on one point Q will remain when integrating 
(or tracing) over the link variables that have one end 
point at site i. 



J DUeM-^UCU]^l[fiC^) 



(84) 



where the 2d{d — 1) fields have been lumped into the 
vector Q. The effective dual action 



Seff = -^Vafl-nsCapCfS - ^ln/(5). 



(85) 



This is a 2d{d — 1) component spin system with pair 
interactions and a nontrivial onsite interaction. As V 
corresponds to nearest neighbor interaction on the even 
sublattice (~ k'^ for small wave numbers), its inverse 
corresponds to Coulombic interactions (~ k~'^) in the 
long wavelength limit. More precisely, the dual model 
corresponds to coupled Coulomb gases having a net of 
2d(d+l) components (in the Abelian case) instead of one. 

The function f{Q) is always even and is further al- 
ways invariant under the group symmetry of U. In the 
argument of the exponent of Eqn.(p4[) 



{9Ug-'gUg-\r 



= Tr 
= Tr 



{gUUg-\r 

{g-'CgUur 



(86) 



for any integer power p > 1 in the expansion, where the 
cyclic property of the trace and the group invariance of 
the measure DU were employed. A pure gauge theory is 
equivalent to gauge invariant charges interacting by pair 
Coulomb potentials. 

As the two dimensional Z2 lattice gauge theory is 
equivalent to the one dimensional Ising chain, its dual, a 
special isotropic two dimensional Coulomb gas with spin- 
spin interactions, has a canonical one dimensional behav- 
ior. (Its partition function Z = [{coshK)^ + (s'mh K)^]^ 
(with L the linear extent of the system) if the two di- 
mensional system has periodic boundary conditions along 
both axis (each of length L)). This is a simple exam- 
ple of a dimensional reduction. The same applies to a 
C/(l) system where in the two dimensional dual model- a 
special, eight component XY type model with Coulomb 
like interactions- we will find that the free energy den- 
sity -N-'^lnZ tends to {-ln[Io{K)]) with Iq a Bessel 
function. 

In fact, by selecting the temporal gauge, in which the 
gauge links Uo{n) = 1 along all links along the time (0) 
direction, any two (or 1+1) dimensional gauge theory 
with plaquette terms can be recast as a classical one di- 
mensional model having nearest neighbor interactions. If 
we dualize any two dimensional gauge theory we obtain 
a special, spatially symmetric, two dimensional Coulomb 
gas that is a nearest neighbor one dimensional system in 
disguise. 

To get an appreciation for what is happening, let us 
examine matters in detail. For the Z2 lattice gauge the- 
ory in c? = 2 dimensions, the vector ^ is four-dimensional 
and at each site. 



f{() — 16 (cos cos ^2 COS ^3 cos C4 
+ sin Ci sin C2 sin C3 sin ^4) , 



(87) 



where we have marked in a counterclockwise fashion the 
bonds {U} about a given site i by Ui,U2,Uz and U^. 
With these bonds we identify 4 dual fields in the following 
way- Ci arises from a coupling to the two bond product 
C/4C/1, C2 arises from a coupling to the product U1U2, and 
so on cyclically. 

The interaction kernel in the 4 dimensional internal 
space of C,(k) (the components of Q in internal and mo- 
mentum space), reads 



v = -K 



A2 
A2 



with the two by two submatrix 

cos(A;i + k2) 



Ao = 







cos(A;i — ^2) 



(88) 



(89) 
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where ki and fc2 denote the x and y components of the 
momentum in the plane. 

As detailed earlier, for all K < Kc, a symmetry oper- 
ation may be performed to map any spin model to onto 
a higher (£>>!) dimensional one, by merely replacing a 
one dimensional momentum coordinate on a spin chain, 
k, by a linear combination X^jLi ^j^j where each of the 
elements {Mj} is a constant matrix, in the argument of 
the interaction kernel. This is exactly what is done here. 
The kernel that appears in our two dimensional problem, 
V, may be obtained from a one dimensional kernel with 
the substitution cos A; — > cos{ki + k2a-3) with cts the Pauli 
matrix to the standard nearest neighbor kernel on a spin 
chain. This transforms the argument of the cosine, k, 
into a two dimensional matrix and makes the interaction 
look two dimensional (involving both ki and fc2) while 
in reality, the system is one dimensional spin chain in 
disguise. 

The interesting question is, of course, whether we 
may see similar dimensional reductions (exact or approx- 
imate) in higher dimensional gauge theories of higher 
groups. 

VIII. PHASE INTERFERENCE AND A 
MAPPING ONTO A SINGLE SPIN PROBLEM 

We will now map classical high dimensional (rf > 1) 
spin problems onto translationally invariant systems in 
one dimension. 

For the nearest neighbor ferromagnet, the basic idea is 
to make a comparison between the the kernels 

d 

Vd = -2^cosfc; 
1=1 

d 

ui = -2^coscifci (90) 

1=1 

in d and in one spatial dinicnisions respectively. Note that 
the one dimensional problem has, in general, d different 
harmonics of the single momentum coordinate fci . If the 
coefficients c/ are integers equal to Ri then the one di- 
mensional problem amounts to a ferromagnetic chain in 
which each spin interacts with d other spins at distances 
{Ri}f^i away. When evaluating the loop integrals, we 
will find that the partition function/ free energy devi- 
ate from their d dimensional values due to "interference" 
between the various harmonics in lii(fci). If all harmon- 
ics acted independently in the effective one dimensional 
problem then the d dimensional result would be repro- 
duced. The advantages of the form Oi are obvious. Per- 
haps the most promising alley is that of "incommensurate 
dimensional reduction". If this method may be also ex- 
tended to fermionic systems that it would suggest that 
we might examine d-dimensional problem described by Vd 



by a bosonization of the one-dimensional problem of vi. 
A solution to the one dimensional problem posed by vi 
with mutually incommensurate coefficients {c^l^L^ will 
immediately lead to the corresponding d— dimensional 
quantities. 

A. Commensurate dimensional reduction 

Consider a one dimensional ferromagnetic lattice 
model with a nearest neighbor (or Range = 1) interac- 
tion and a Range = n interaction: v{k\) = — 2(cosfci -|- 
cosnfci). The loop integrals containing terms of the 
cos fci origin and terms stemming from cos nki are in- 
dependent to low orders. Mixed terms can survive only 
to orders (l/T)""*"^ and higher. The delta function con- 
straints X^Qi = would generate terms identical to 
those of the two dimensional nearest neighbor model with 
v(ki, k2) = — 2(cosA:i -|- cosfc2) with the two delta fimc- 
tion constraints for the two separate components of the 
momentum, J^Qi = ^ and J^l'^ — 0' at every vertex. 
Thus, in a sense, this system is two dimensional up to 
T"". Similarly, if a one dimensional chain has interac- 
tions of length one lattice constant, R = n and R = m 
i.e. 

v{ki) = — 2(cosfci -|-cosnfci -|-cosmfci) (91) 

then this system, as fleshed out by a 1/T expansion for 
the partition function or for the free energy, is three di- 
mensional up to order = n if m = n^. (If m = n -|- 1 
then the triangular ferromagnet is generated). Similar, 
higher dimensional extensions similarly follow from the 
lack of commensurability of the cosine arguments (for 
four dimensions 

v{ki) — — 2(cosfci -|- cosnki + cosmfci -|- cossfci) (92) 

with s = n^ etc.) 

The corresponding transfer matrices are trivial to write 
down. For the three dimensional case: 

{si...Sn2\T\s[...s'^2) = exp[/?{i(siS2 + ••• + S„2_iS„2) 

+ + - + Sn2_i<2) 

+ 2(*1*"+1 + S2Sn+2 + ■■■ + S„2_„S„2) 

+ 2 ~^ ^2-^n+2 + ■•■ + ^n^-n.-^n'^) 

+ Sn2s[ + {Sn2-n+is'i + S„2_„_,_2S2 + -■ + 3^2 s' J 

+isis[ + ... + s„2s;2)}] 

= eMf3'Tji]. (93) 

In higher dimensional generalizations, Tji will be slightly 
more nested with the span replaced by n'^~^. In 
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Eqn. (|93|) , the matrix indices i and j are written in bi- 
nary numerals in terms of the spins 



(94) 



where a is the one dimensional coordinate along the bra 
spin indices. By Eqn.(p3|), Tji are expressed in terms of 
a sum over two digit products where the digits are those 
that appear in the binary (length n spin) representation 
of the coordinates j and i. 

As usual, the partition function 



Z = Tr[T^] = E Af 



(95) 



where {A^} are the eigenvalues of the transfer matrix. In 
the following A„iax will denote the largest eigenvalue. 

Note that the transfer matrix eigenvalues correspond 
to periodic boundary conditions as strictly required: we 
employed translational invariance to write the partition 
function expansion in Fourier space. If the boundary 
conditions are not periodic then all of this is void. 

By the lack of interference effects, the connected dia- 
grams for the free energy are correct to 0(/3") for inter- 
actions with m — (Eqn.^l]). This along with 



0f = J2fp^' 



InZ 



In An 



(96) 



(where the last equality holds in the limit of large sys- 
tem size N) and a canonical power expansion having an 
infinite radius of convergence 

oo 

A„.ax = exph E fpPP] = E (;^(- E fp^'^y 

p m—0 ' p 

= J2b,(3P (97) 

p 

imply that A^ax is correct to O (/?"). 

If we already know the lower order coefficients corre- 
sponding to 0{l3P) with p = 0, 1, 2, 3, {n — 1) when we 
consider the Range = n problem, then we may set out 
to compute is the coefficient of This will amount to 
extracting only a single unknown (i.e. the term corre- 
sponding to the coefficient of /3") from a transfer matrix 
eigenvalue equation, This can be done to the next order 
etc. recursively. A simple yet very long single linear rela- 
tion gives the largest transfer matrix eigenvalue to each 
higher order in the inverse temperature p. 

Explicitly, we can write longhand 



det(T- A) 



nK-.,(/3)-A(/3)<5,,J = (98) 



All one has to do is to yank out the 0(/3") term from 

2 

the sum over products of 2" terms. If we pull out a 
power pji- from each of the elements [Tji. {(3) — X{P)Sji.] 
then those will need to satisfy 



(99) 



In each element of the transfer matrix, exp[/37^i], the 
coefficient of is trivially Tj^/pl. For each permutation 
(or "path" of coordinates (i,j) within the matrix) to be 
summed for the evaluation of the determinant 



e 

path 



n (100) 



is exactly the product of the coefficients of /J^'^ such that 
those lead a net power /3" (i.e. satisfying Eqn.(|99|)) from 
that particular path. 

It follows that net contributions stemming from Tji . (/3) 
in Eqn.(||) are 

^^"^^^ E ^p-'^t E (101) 

all paths U'^j) given path 

where Cpath simply denotes the sign of the given permu- 

tation ^ {ij}. 

To take into account the products including X{/3) let 
us define 



D 



(n:p) 
2m 



1 



(n-p) 



E 



all 2m paths 



^path[ -^J 
(jij) in a given path 



n—p 



(102) 



where the summation is over paths going through all 
possible {2m) given points on the diagonal and de- 
notes a summation over Tji^ in the given paths sans the 
contributions from the (2m) diagonal points. The de- 
terminant paths in J^aii 2m paths '^^^ p&ss through more 
then (2m) diagonal points- it is just that we need to sum 
over all those that in their pass also traverse all possible 
sets of (2m) given points on the diagonal and where those 
points are excluded from the second Tji^ summation. 

The /?" component of Eqn.(p8|) is 

n 2" "1 

Fi-)+Y^J2 Dt"^ E V-^P2™=0 (103) 

where {^p}p=o appear in the expansion of the largest 
eigenvalue (Eqn.(p7|)). 

If the coefficients {bpYpZ^ then Eqn.( p^ ) where 6„ ap- 
pears there only linearly gives 6„ and one may then pro- 
ceed to look at the eigenvalue equation for the next value 
of n. 
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In the absence of an external field, for n = 1 the 
largest eigenvalue for the d dimensional problem is A+ — 
2cosh(i/3 and other eigenvalue A_ = 2sinhd/3. For 
the three dimensional case setting m = = n = 1 
leads to the one dimensional momentum space kernel 
v{ki) = — 6cos/ci which in real space corresponds to 



H — —3 SiSi+i. 



(104) 



Of the two eigenvalues only A-|_ has a nonzero bo- Know- 
ing the values bo = 2 and 6i = of the l arge st eigenvalue 
we may proceed to find 62 from the Eqn.(103) when n = 2 
etc. This, of course, may be extended to systems with 
magnetic fields. 



Setting bo = 2 in Eqn.(103) wc find 



2" " 

bn J2 ^^""D. 



m — 1 



An) 



71-1 



+ J2 J2 bp,bp,D[ 



(ti;p) 



P=0 p=pi+P2 



n 2" -1 
p=0 m=2 



E 

^^."^^Pi— P with Piy^P 



(105) 



This is the rather explicit recursive relation giving 6„ 
once {bp}pZo are known. To find out the coefHcient of 6„ 
let us look at 



exactly 2s diagonal paths 



(106) 



where the summation admits only those permutations 
that lead to exactly 2s diagonal elements. For the 2" 
diagonal path in the matrix T: epath = 1- For all paths 
threading 2"^"^ diagonal sites Cpath = — 1 and the num- 

ber of such paths is (^^ 
quantity amounts to 



^path 

For the general case the last 



Evaluating, we find 

2" 



{-ly 



2" 
2s 



(107) 



E D'r' 



(2" 



r r 

'^Fpq[{l, — 



(4m)! (2" - Am)\ 
1 m 



Am ^ , 

{3 + Y'l + 2™},1] 



(108) 



where Fpq is a generalized hypergeometric function. 

Insofar as simple geometric visualization is concerned, 
it is amusing to note, as shown in Fig.(P) for the two 



3n+l 3n+2 3n+3 
• • • 



2n+l 2n+2 



4n 




FIG. 3. Lacing of an (n x L) slab (with L 00). All 
real space high temperature diagrams up to order (tanh" (5) 
such as the closed rectangle shown above are of equal value 
in both this one dimensional system with interaction ranges 
i? = 1, n and the full (L x L) two dimensional model with near- 
est neighbor interactions. The width of the slab "n" serves as 
an inverse temperature axis in the following sense- the wider 
it is, the higher order in l/T that we may advance towards 
the full two dimensional model. We may get an analogous 
"three-dimensional" slab with "two imaginary time axis" if 
we set m — . The path will then thread sites before 
continuing upward. 



dimensional case, that employing the usual high temper- 
ature expansion in powers of tanh /? one would reach the 
same conclusion regarding the correctness (up to order 
n) of the partition function evaluated for a width n slab 
vis a vis the partition function of the two dimensional 
system. We may look at the a finite thickness (n) slab of 
the two dimensional lattice along which we apply peri- 
odic boundary conditions. Let us now draw a string going 
along one row of length n, after which it would jump to 
the next row, scan it for n sites, jump to the next one, 
and so on. On the one dimensional laced string, the sys- 
tem is translationally invariant and the interactions are of 
ranges R— l,n. Or, explicitly, by counting the number 
of closed loops in real space (employing the standard, 
slightly different, diagrammatic expansion in powers of 
tanh /?), we see that the terms in this expansion are also 
identical up to order tanh" /3 |^l|] . 

A lacing of a two dimensional slab by a one dimensional 
string, as shown in the figure, is one of the backbones 
of Density Matrix Renormalization Group Theory when 
applied to two dimensional problems. 



B. Incommensurate dimensional reduction 



The one dimensional continuum limit (—A < fci < A 
with A — > 00) momentum space kernel 
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cos ciki 



(109) 



1=1 



with mutually incommensurate {c/} will reproduce the 
expansion with the kernel in Eqn.(^ for the d dimen- 
sional hypercubic lattice system ||22 1 . 

The proof is quite easy. The most general integral is 
of the form 



n dkb 

A .-A 



b=l 



propagators d loops 

^exp[z ^ Mabhci] 

a=l 1=1 b=l 



(110) 



where the index b runs over the various independent loop 
momenta (k is still merely a scalar for this one dimen- 
sional problem). Unless, for a given fcf, integration, the 
argument of the exponent is identically zero (i.e. corre- 
sponding to a term that would be generated in the d- 
dimensional nearest neighbor problem) an "interference 
term" results. However, such a term is down by ^^(A"^) 
by comparison to the "good" noninterference terms that 
occur in the d-dimensional problem. The canonical loop 
integral reads 



1 

2A 



dkh exp[zc/ch 



TrS{c) 
A 



(111) 



where c is a linear combination of the d coefficients {q}. 
Unless c vanishes identically (i.e. unless a "good" non- 
interference term occurs), then c ^ if {q} are chosen to 
be mutually incommensurate. Whenever c = 0, the in- 



tegral in Eqn.(lll) transforms into the Kronecker delta 
(5(c, 0) = 1 in the limit A — > oo and the integral at- 
tains exactly the same value that it would take on for 
the d— dimensional hypercubic lattice system. Similarly, 
in this limit all "bad" interference terms will evaporate 
in all diagrams. 

As this holds for any diagram, this conclusion is also 
valid for any set of diagrams (not only those encountered 
in computation of the partition function or the free en- 
ergy), e.g. Cy, X etc. 

As a consequence we may state that for a single spin 
1/2 particle with an action 



dio 'ip{—uj) v{uj) 1p{Lij) 



(112) 



where i}}{t) is the two component spinor, the partition 
function is identically the same as that of the d dimen- 
sional nearest neighbor ferromagnet. 

The proof of this statement trivially follows from 
breaking up the [0, (3] segment along the imaginary time 
axis into L pieces and allowing L ^ oo. 



We have mapped the entire three dimensional Ising 
model onto a single spin 1/2 quantum particle p^ . 
The Fourier transform of Eqn.(10£) on a lattice reads 



V{xi) = -2^[(q -t-a;i)sin7r(c; ~ xi) 
1=1 

+ (c/ - xi)] sin7r(a::i + q)] 
x[(Q-Xi)(Q+a;i)]-i, (113) 

a sum of shifted Coulomb like interactions. If an aver- 
age (2A)^^ dci is performed over each of the d coef- 
ficients {c/} then the resulting quantities are the exact 
d-dimensional ones. If "bad" interference terms result 
then their value is 0(l/c). In the averages, it is the large 
coefficients that dominate and these will lead to a van- 
ishing average value of all interference terms. 



IX. INCOMMENSURATE DIMENSIONAL 
REDUCTIONS OF QUANTUM MODELS 

The incommensurate dimensional reduction averaging 
theorem which we just proved by interchanging the order 
of integrations over {q} and the loop momenta {kb\ holds 
for any translationally invariant system (spin, bosonic, 
fermionic) in which the unperturbed action 5*0 links two 
nearest neighbor sites. For a {d+1) dimensional quantum 
lattice system, the Green's function conjugate to Sq in 
the undualized (non Hubbard-Stratonovich transformed) 
system 



Go = 



.( 



d+l 

2 cos fc; + A 
1=1 



(114) 



As before, if the magnitude of the squared mass \A\ = 
const > 2(c?-|-l) then a geometric series expansion in A^^ 
may be performed to generate a sum of harmonics. For 
each given single diagrammatic term, our previous proof 
may be applied. 

We note that fermionic electron (or other) fields might, 
perhaps, be formally bosonized on each individual chain. 
The one dimensional band, 



Vl 



-2^ cos ciki 
1=1 



(115) 



is now a simple sum of d cosines (au lieu of the standard 
single tight binding cosine). For each given set {c/jf^j^, 
the band dispersion of Eqn.(115) may be easily linearized 
about its respective Fermi points. Consequently the stan- 
dard bosonization methodology may be applied. In the 
large {q} limit, the number of Fermi points becomes in- 
finite. A one dimensional Fermi system with (2n) Fermi 
points is expected to be equivalent to a d-dimensional 
system up to 0(/3"). 
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The astute reader will note that we may just as well 
use Eqn. ( 114 ) to study a classical spin system (with the 
upper bound d + 1 replaced by d) and prove everything 
that we have previously stated without the need to con- 
sider the dual model. Such an approach is not as rigorous 
as that followed hitherto. First, the geometric series ex- 
pansion in the harmonics exp[ifc;] converges only for a 
sufficiently absolute value of the mass- or, equivalently 
for temperatures far enough from the mean field temper- 
ature. Second, the expansion in Go is in general explicitly 
asymptotic as the perturbing piece is of higher order in 
the fields than Hq. The dual formulation was free of these 
pathologies. 



X. PERMUTATIONAL SYMMETRY 

The classical spin spherical model (or 0(n oo)) par- 
tition function 



and consequently 



Z — const 



1 



(116) 



where the chemical potential /i satisfies 

/3: ' ' 



is invariant under permutations of {v{k)} 
the above, the permutations 



N 



{Pk} 



(117) 



{v{Pk)}. In 



(118) 



correspond to all possible shufflings of the N wavevectors 
ki. Although quite simple, this is not universally realized. 
Several authors have attempted to compute the critical 
exponents in the spherical limit (via an RG calculation) 
for systems having different minimizing manifolds yet all 
sharing the same relevant density of states. This quest 
was not very economical. As unrealized by these authors, 
by permutational symmetry these models are identical. 

This simple invariance allows all d-dimensional trans- 
lationally invariant systems to be mapped onto a 1- 
dimcnsional one. Let us design an effective one dimen- 
sional kernel Ve//(fc) by 



S[v{k) - v]d'^k 



(119) 



The last relation secures that the density of states and 
consequently the partition function is preserved. For the 
two-dimensional nearest-neighbor ferromagnet: 



.dV, 



eff 1-1 



dk 

dx 



Cl 



'o VT^^l-{Veff+X-2)^' 



(120) 



HVeff) 



Cl 



F(sin" 



(3 — u)u ' 



-)du, 



(121) 



where F(t, s) is an incomplete elliptic integral of the 
first kind. Eqn. (121) may be inverted and Fourier trans- 
formed to find the effective one dimensional real space 
kernel Veff{x). We have just mapped the two dimen- 
sional nearest neighbor ferromagnet onto a one dimen- 
sional system. In a similar fashion, within the spherical 
(or equivalently the 0{n — > oo)) limit all high dimen- 
sional problems may be mapped onto a translationally 
invariant one dimensional problem. It follows that the, 
large n, critical exponents of the d dimensional nearest 
neighbor ferromagnet are the same as those of transla- 
tionally invariant one dimensional system with longer 
range interactions. We have just shown that a two di- 
mensional 0(n 3> 1) system may has the same thermo- 
dynamics as a one dimensional system. By permutational 
symmetry, such a mapping may be performed for all sys- 
tems irrespective of the dimensionality of the lattice or 
of the nature of the interaction (so long as it translation- 
ally invariant). This demonstrates once again that the 
notion of universality (with dependence only on the or- 
der parameter symmetry, dimensionality etc.) may apply 
only to the canonical interactions. 

The lowest order term breaking permutational sym- 
metry in our high temperature expansion is rf{x)rf{y). 
Thus permutational symmetry is broken to 0{/3^) for fi- 
nite n. For a constraining term (e.g. X^x 
for 0(1) spins) symmetric in {ri{k)} to a given or- 
der, one may re-arrange the non-constraining term 



J2i:v-\Pk)\7]{Pk)\^) and rela- 



bel the dummy integration variables H[{r](k)}] — > 
H[{r]{P~^k)}] to effect the constraining term augmented 
to a shuffied spectra v{Pk). [ p4[ 
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XI. APPENDIX A: THE HUBBARD 
STRATONOVICH TRANSFORMATION FOR 
0{N) SPIN SYSTEMS 

For the benefit of our uninitiated readers we present 
the standard Hubbard Stratonovich transformation. 
This "transformation" merely relies on the elementary 
Gaussian identity 
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-^cxph^AV + 077] = A-i/2exp[^]. (122) 

As well known, this easily generalizes (upon diagonaliza- 
tion) to 

J {Dt]) exp[-i(?7, Kri) + (77, 4>)] 

= {detK)-^'^ eKY,]^{(t),K-^(t>). (123) 

For our particular purposes, note that 

v/det(2^(-y-i)) = j jj^ exp[-i(sy-77) 

{-V-'){Vs~7j)]. (124) 



Thus 



Z = ^exp(-is,ys) 



det (— ^) y exp[i77, V S] 
^exp[-s,77] 



det(— / L'7/exp[-iJ] 



where 



H = -J2 ln(cosh7;,) - i(7,, V-'v). 



(125) 



(126) 



In the above s = (5*1, S2, ■■■tSn) where i labe l the sites 
X. For an Ising system S{x) — ±1. In Eqn.(126), the 
relation 

exp(-s, 7;) = ^ exp(-si77i - ... - snVn) 

AT 



]^(2cosh7/,) (127) 



z=l 



was employed. 

In an 0{n) spin system, tracing out over the spins leads 

to 

^ exp(-s, ??) = H [ / dns, exp(-S^i • fft) 

is) »=1 

^ r 2 1 

= n r(7i/2) {—Ti^-^ ini2-A\m)\) ■ (128) 



XII. APPENDIX B: THE EIGHT-FOLD 
"PORCUPINE" SPIN MODEL AND OTHER 
"PLATONIC" SPIN MODELS 

For "Z2(8)-^2® -2^2" models, the three-dimensional spins 
explicitly read 



with 



S = (CTI, (72, (T3) 



= ±1 



(129) 



(130) 



(i.e. spin \ like). This corresponds to the 8 possible 
states (the vertices of the unit cube) (±1,±1,±1), or in 
a more pompous notation- the eight direct product states 
of three spin-1/2 particles \a\) ® \a-i) (g) jcra) states. For 
most canonical rotationally and translationally invari- 
ant models {(y\aio^\R\a\aiOj) decomposes into a sum of 
three decoupled Hamiltonians each living within its own 
individual two dimensional (|cri)) subspace. The three 
dimensional "single particle (spin)" problem will reduce 
to a sum of the three individual non-interacting Hamilto- 
nians each corresponding to a single free (spin-1/2) "par- 
ticle" having only two states. A Hamiltonian containing 
only [Si • Sj] terms (with i and j lattice sites), amounts 
to 

3 

i?({Si})= E nxi,Xj)[Si.Sj] = ^7^({af}) 

Xi,Xj Q— 1 

(131) 

where 7i are Ising Hamiltonians, i.e. 



H(K"})= E < nxi,X2)a2° 



(132) 



The partition function 



z=l[z^ = {Zi 



(133) 



Q = l 



The Z2 model behaves exactly like an Ising model (or a 
q = 2 Potts model). This may be extended to other N 
state clock models. When N = 2 the "clock model" is 
trivially the q = 2 Potts (Ising) model. When iV = 77. -I- 1 
in the 7i— dimensional "clock model", a mapping to an 
(77,-1-1) state Potts model is possible (by constructing 
a regular 77, dimensional tetrahedra). For N = 2" (the 
77 dimensional hypercube), the model may be mapped 
(once again) onto 77, decoupled q — 2 Potts models. In 
two dimensions, = 4 is the critical number of states 
required for the nature of the transition to vary in the 
unfrustrated ferromagnetic Potts model. Thus, in the 
unfrustrated N is monotonically increased (for 

77 > 3), the nature of the transition changes, at least 
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twice, from that with {N = 2; 2 = q < Qc) to {N = q = 
n + 1 > (7c) to the q < qc case {N — 2"; 2 = q< q^) once 
again. For all < n + 1 uniformly spaced polarizations, 
with overall isotropic orientation, the model is exactly 
equivalent to & q = N state Potts model; in this case, the 
discrete polarizations span an A^— dimensional subspace 
embedded in the full n— dimensional spin space. 

The three-dimensional "clock model" is restricted only 
to very special N values corresponding to the perfect 
solids (the six perfect polyhedra - the "Platonic solids" - 
when n = 3) and when n > 3, the uniformly spaced spin 
polarizations are restricted to the vertices of the (special 
small set of allowed) regular polytopes possible for that 
value of n. The N ^ oo limit may not be taken- non- 
infinite iV's are bounded by a finite number (e.g. iVmax = 
60 for 71 = 3). 

When N = oo (not the limit) for a given n the con- 
tinuous 0{n) model is attained -this model has second 
order phase transitions in d = 3 for its unfrustrated fer- 
romagnetic variant. 

When n — 3, the only possible finite N variants are 
A'^ = 4 (the tetrahedron spin model that may be mapped 
onto the q = A state Potts model as just discussed), 
N — 8 (the cube or hexahedron - gg// = 2 as discussed 
in on page 1), N — 6 (the octahedron), iV = 20 (the 
dodecahedron), iV = 12 (the icosahedron) , and TV = 60 
(the soccer-ball or truncated icosahedron). 

For the iV = 6, 12, 20, and 60 variants a mapping into 
a model similar to the Potts model is possible. Here the 
number of allowed individual two-spin energy states is no 
longer effectively two (as in the Potts models for equal 
and unequal spin polarizations of neighboring spins) but 
is slightly larger - a trivially generalized Potts model (e.g. 
there are three different possible values for a spin-spin 
interaction in the iV = 6 model). 

The octahedron {N = 6) model is equivalent to a 
"three-color dilute Ising model" where each site is col- 
ored red, white, or blue. On each individual color clus- 
ter independent Ising spins interact amongst themselves. 
Consequently, a sum over all possible colorings (dilute 
inter penetrating color sites) of the lattice sites is to be 
performed. 
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V. V. Bazhanov. R. J. Baxter. J. Stat. Phys., 69 453 



(1992), |hep-ph/9212050| 
Y. Tokura and N. Nagaosa, Science 288, 462 (2000), and 
references therein 

J. W. Negele and H. Orland, "Quantum Many-Particle 
Systems" Frontiers in Physics, Addison- Wesley (1988) 
R. J. Baxter "Exactly solved models in statistical 
physics". Academic Press (1982) 

As a technical aside, note that on the lattice (unlike 
the standard continuum theory) each vertex conserves 
momentum only up to a reciprocal lattice vector K = 
27r(ni, n2, rid) with rii integers. Formally, this is as x 
only take on discrete values on a lattice. 



^ exp[ifc • f] = iV ^ 



(135) 
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where K are reciprocal lattice vectors. 
Thus at each vertex of order 2m we have 



(2m)! ^ 

K 



(136) 



where qi„ denotes the net momenta flowing into the ver- 
tex. As for all lattice models the propagators are periodic 
v{q) = v{q + K) for each of the N reciprocal lattice vec- 
tors K. As such the momentum conservation constraint 
may be replaced by 



1 



(2m) 



K,Q)^Nj^&{q,^,Q). (137) 



In the aftermath, each vertex carries a factor of the vol- 
ume (A'^) multiplying the discrete Kronecker delta func- 
tion. After the momentum conservation constraints are 
taken care of at each vertex, we will be left with a summa- 
tion over all remaining independent loop momenta {fct,}. 
Each one of these summations is to be replaced by the 
standard 



E 



N 



(Tk 
{2-kY 



Thus, at the end we will be left with a factor of 
^loops -I- vertices - propagator lines 



(138) 



(139) 



multiplying an integral over the independent loop mo- 
menta {fci,} in Eqn.(p^). Compounding all together, as 
in each closed diagram the number of independent loops 
+ number of delta function momentum conserving ver- 
tices = 1-1- number of propagator lines, a net factor of A'^ 
should be attributed to each independent closed bubble 
just as it does in the conventional continuum theories. 
[16] As noted earlier, for any lattice model 



v{k) =v{k + K), 
as exp[iK ■ x] — 1 in 



v{k) 



V{x) exp[ifc ■ 



(140) 



(141) 



The sum over x spans all discrete lattice sites. 
The way to analytically continue v{k) when fci becomes 
complex is to make the k^ dependence also invariant un- 
der shifts by {2n). Within each Brillouin Zone the 
dependence stems from the analytic continuation of the 
functional dependence on the real part. The only compli- 
cation is the Gibbs' phenomenon in kg which occurs when 
taking the Fourier transform of the lattice V{x). This is 
slightly avoided for the real part of V{x) as follows: Un- 
der parity and all individual reflections: V{xi,X2,xs) — 
V{xi,X2, —xs) etc. This implies that in the Fourier trans- 
form v{k): v{k\,k2,kz) — v{ki, k^, —kg). Thus, taking 
only the even part in fea of cos(fc2 + iXks) we will only 
be left with the terms cos k2 cosh fes. Now, cosh fea is con- 
tinuous and periodic when the Brillouin Zone is chosen 
as —TV < ka < n. No Gibbs phenomenon occurs. The left 
and right derivatives are still discontinuous at the zone 
edges. 



[17] By replacing the scalar product of Eqn.(|l|) by a non- 
rotationally invariant product of the form 



5(f)oS(y)=^[5.(5)]^[S.(y)]^ 



(142) 



with a power z 1, the rotationally invariant Bessel 
functions (with prefactors) of Eqn.(^ are replaced by 



n p 

F{fi) = ^ / dCl rja exp[iriaCll 



(143) 



with Clct denoting the a component of the n-dimensional 
unit vector Q. The rotational invariance of the Bessel 
function terms in Eqn.(H) is replaced by the point group 
symmetry of F{'if) (that of the n-dimensional hypercube). 
By repeating the steps leading to Eqn.([70|), we find that 
the partition function of 



H = 



- E ^'^[^' 

(ij) along (1,1) 

E ^'^[^' 

(ij) along (1,-1) 



(1)1 



(1)1 



(2)iz [^(2)12 



(144) 



is identical to that of a one dimensional system with 

= oSj. (145) 



[18] P. Fazekas, :ond-mat/0011354 



[19] E. Berg, E . Altman, and A. Auerbach, cond 



mat/0206384 



[20] For the (7(1) case, the real action 



Sr = -KRe{^{UUUU)}. 



(146) 



As usual, in the non-Abelian case, the real part of trace 
of the plaquette product must be taken. 
[21] It might be easier to compute the terms in the (3 ex- 
pansion for a three dimensional Ising nearest neighbor 
ferromagnet by going to a triangular lattice and delet- 
ing all bad "constructive interference" loop integrals. A 
line threading the triangular lattice will connect a spin 
with spins of distance Range = 1, n, and n + 1 along 
the line. "Bad interference terms" can result from loop 
integrals of the form {exp[ik] exp[mA:] exp[— i(n + l)fc])*'. 
This is the only way in which loop integrals can conspire 
to give nonzero results when n — + oo. Thus, trivially, the 
high temperature expansion for d = 3 cubic ferromagnet 
is equivalent to that of the triangular ferromagnet up to 
order /3^~^. Similarly, if we set s = n — 1 we obtain a 
2-Dimensional lattice model in which each spin interacts 
with its four nearest neighbors and four diagonal neigh- 
bors. 

[22] In the continuum limit, this corresponds to the kernel 



or, in a lattice approximation (— tt < fci < vr), to 



(147) 
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d 

V{xi) = -2^[(ci + a::i) sin7r(c; - Xi) 
1=1 

+ {ci - xi)] sin7r(a;i + c;)] 
x[{ci - xi){ci + xi)]-\ (148) 

(i.e. a scenario in which each spin is effectively com- 
posed of d shifted "Coulombic sources" (more precisely, 
each spin is composed of d pairs of "charges" generating 
shifted sine potentials)). 
[23] The astute reader will note that the partition function 
of the d— dimensional nearest neighbor Ising model on a 
lattice (with ~n < ki < tv) is equivalent to that of the 
d— dimensional nearest neighbor Ising model in the con- 
tinuum (with the same kernel v{k) but with unrestricted 
(-A < fe < A)) 

[24] That permutational symmetry cannot, for finite n, hold 
to all orders in /3 immediately follows from the fact that 
the normalization constraints of the 0{n) spins are not 
permutationally symmetric. If 

then the normalization constraint at each site x reads: 

1 = ^^(f) = -i- V s(fe) ■ s{k') exp[j(fc -f- fe') • x\. (150) 
jV^ z — i 

k,k' 

A permutation of {s*(fe)} (which is equivalent to a per- 
mutation of the spectra {v{k)} within the quadratic in- 
teraction kernel) would transform the A*' constrains of 



Eqn.(15[) into another set of equations. 
The sole normalization constraint which is permutation- 
ally symmetric is the global normalization constraint (ap- 
propriate to the 0(n oo) spherical model) 

l = ^$]5"'(fc)-s(fe"*'). (151) 

fe.fc' 

This illustrates why permutational symmetry must triv- 
ially hold, to all orders, in the 0{n oo) limit (and only 
in this limit). 

For the quantum spin system, 

s{k)xs{k')=iNsg^p. (152) 

Just as with the normalization conditions for the clas- 
sical spin models, these commutation relations are not 
invariant under permutations of {s(fe)}. 
We mention, in passing, that permutational symmetry 
trivially holds in all quadratic (non-interacting) transla- 
tionally invariant Bose and Fermi systems. The respective 
commutation and anticommutation relations 

[at(fc'), a(fc)] = {c\k), c{k')} = 4.fe' (153) 

are invariant under permutations 0{k) — * 0{Pk) for all 
bosonic or fermionic O. 
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